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SYMBOLS 


A 

b 


c^,  c2 


E[U  ] 
E,  E 


f,  fx,  f2 


amplitude;  mean  amplitude;  inner  radius  of 
hollow  cylinder 

amplitude 

outer  radius  of  cylinder 
outer  radius  of  case 
arbitrary  constants 

exponential  function 
variance  of  temperature 

modulus  of  elasticity  of  cylinder,  of  case,  psi 

average  frequency,  cy/hr 
frequencies,  cy/hr 


fAU) 

f 

m 

8m 

h(.) 

h 

c 

i 

j 

V')’  J<‘>01 

k(-) 

l(-) 

m 


density  function  of  amplitudes 
Kelvin  function 

Kelvin  function 

function  h(‘)  * ber(-)  + i bei(*) 
thickness  of  case 

index  (j  = r,  9,  z) 

Bessel  function  of  the  first  king,  order  zero; 
real  and  imaginary  parts 

function,  k = ber(-)  - bei(') 

function,  l = ber(*)  + bei(*) 

order  of  Kelvin  function 
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p 


R(.) 


S(.), 


S 

m 


t 

T 

U(r ) , U , U 
m y 


W 


i’ 


W„ 


function  in  stress  equations 

Probability 

radial  coordinate,  in. 

frequency  response  function  of  temperature 

power  spectral  density  functions;  for  input  and 
output  temperatures 

radial  transverse  and  axial  stresses,  psi 

mean  stress,  psi 

time,  hrs 
period,  hrs 

temperature;  annual  mean;  amplitude  of  annual 
cycle 

one  sided  power  spectral  density  function; 

2 

for  input  and  output,  °F  /cy/hr 


x 


variable 


x = 


r 


Y Y Y 
O’  OR’  01 

z 

a 

a,  ac 

r 

S 

♦ 


n 


Bessel  function  of  the  second  kind  of  order 
zero;  real  and  imaginary  parts 

axial  coordinate,  in. 

2 

thermal  diffusivity,  in. /hr 

thermal  coefficient  of  expansion  in  cylinder; 
in  case,  in./in./°F 

gamma  function 

differential  operator 

phase  angle 

Poisson's  ratio  for  cylinder;  case 
3.14 
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INTRODUCTION 


r 


Engineering  structures  are  subject  to  environmental  variations 
of  the  ambient  temperature  that  in  turn  produce  thermal  stresses  and 
strains. 

Temperature  variations  consist  of  cyclic  components  (annual  and 
diurnal  cycles)  and  random  components  produced  by  clouds,  winds,  rain, 
and  snow.  Solar  radiation  creates  additional  temperature  changes. 

The  response  of  a large  structure  or  structural  component,  having 
various  thermal  properties,  will  depend  on  the  frequency  of  such  changes. 
The  temperature  at  the  interior  of  the  structure  will  be  delayed  and 
attenuated  compared  to  surface  temperatures. 

A slow  variation  such  as  the  annual  summer-winter  cycle  penetrates 
throughout  the  structure  while  the  diurnal  cycle  effects  only  surface 
layers  or  thin  structural  components.  As  a result  temperatures  in 
structures  are  usually  nonuniform  and  produce  nonuniform  cyclic  and 
random  thermal  stresses. 

In  the  following  a large  diameter  circular  cylinder  with  casing 
will  be  examined.  The  annual  thermal  cycle  will  be  considered  to  be  a 
slowly  varying  mean  while  the  diurnal  cycle  and  the  superimposed  random 
variations  are  assumed  to  be  narrow  band  excitations  with  a central 
frequency  of  one  per  day.  The  large  difference  between  the  frequencies 
of  the  annual  and  diurnal  cycles  makes  it  not  unreasonable  to  consider 
the  process  to  be  stationary. 

The  mean  and  standard  deviation  of  the  thermal  stress  at  any  point 
in  the  cylinder  will  be  calculated  from  the  given  mean  and  standard 
deviation  of  surface  temperatures  which  are  assumed  to  be  Gaussian 
distributed.  Because  the  governing  differential  equations  for  tempera- 
ture conduction  and  for  thermal  stress  are  linear,  the  distribution  of 
thermal  stress  at  a point  is  also  going  to  be  Gaussian. 

Once  the  parameters  of  the  distribution  are  calculated  a safety 
factor  analysis  may  be  performed  if  the  distribution  of  the  material 
strength  is  known.  Because  loads  are  repeated  the  safety  analysis  will 
also  lead  to  a life  prediction  capability. 

II.  THE  FREQUENCY  RESPONSE  FUNCTION  OF  TEMPERATURE 

The  frequency  response  function  is  defined  as  the  response  to 
a sinusoidal  input  of  unit  amplitude.  For  the  case  at  hand  the  input 
is  a cyclic  temperature  variation  at  the  surface  of  the  cylinder  while 
the  response  may  be  either  the  resulting  temperature  at  an  interior 
point  or  any  of  the  thermal  stress  or  strain  components  at  the  point. 


Preceding  page  blank 


7 


A long  axisymmetric  cylinder  will  be  analyzed  for  which  the  one 
dimensional  heat  conduction  equation  in  cylindrical  coordinates  is 
valid  [1] 


52u  . I Su  = 1 du 
^ 2 r 5r  a St 


(II-l) 


where  U(t,r)  is  the  time  and  space  (radial)  dependent  temperature  and  a is 
the  thermal  diffusivity  in  in.2/hr  (m2/S). 

The  differential  equation.  Equation  (II-l),  is  subject  to  the 
following  boundary  conditions 

BC.l  : U (b,t)  = eiut 
BC.2  : U (0,t)  = finite 

That  is  a sinusoidal  temperature  input  with  circular  frequency,  o 
radians/hr,  and  unit  amplitude  is  applied  at  the  outer  radius,  b,  of  a 
solid  cylinder  [2], 

The  solution  to  Equation  (II-l)  using  separation  of  variables  may 
be  written  as 

U(r,t)  = R(r ,<d)  eijt  (II-2) 

where  R(r,'j)  is  the  frequency  response  function  of  the  temperature. 
Substituting  into  Equation  (II-l) 

R" (r ,'o)  + ~ R'  (r)  - ^ R(r,w)  = 0 (II-3) 


is  obtained.  Equation  (II-3)  is  a Bessel  differential  equation. 


•>» 

The  above  problem  has  been  solved  by  Dahl  as  far  back  as  1924  for 
a caseless  cylinder  with  the  use  of  Bessel  and  Neuman  functions [3] , 

Here  Kelvin  functions,  which  are  essentially  Bessel  functions  with  com- 
plex arguments  , will  be  used.  Equation  (II-3)  is  satisfied  by 


R(r,w)  = C1JQ 


Vo 


(II-4) 


These  functions  are  described  by  Wattson  [4]  and  by  the  Handbook  of 
Mathematical  Functions  [5], 


In  series  form  Bessel  functions  of  the  first  and  second  kind  are 
written  as 
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Jo(x) 


l 

n=l 


(±M 

(n ! y 


2n  I 


(II-5) 


Yq(x)  = Jq(x)  In  x 


■I 

n=l 


(-0"  M 

(n!)2 


2n 


u 

X 5 


n=l 


(II-6) 


and,  in  this  case,  consist  of  real  and  imaginary  series  Substituting 
N/-i  = ^/i~*  ~ “ into  Equations  (II-5)  and  (II-6)  and  using  x = 7^"r  as 

Jl  a 

the  variable  the  real  and  imaginary  parts  become  the  Kelvin  functions 

2 

(!) 


■V*>  -£■  (-l>" 


and 


2n 


(2n)  ' 


= ber  x 


(II-7) 


(!) 


(2nfl) 


= bei  x 


(II-8) 


(2n+l)! 

The  Bessel  function  of  the  second  kind  may  also  be  developed  analogously. 


It  should  be  recognized  that  y = In  i and  hence 


ln  r)  = 


t1  + in  r 


(II-9) 


Substituting  into  Equation  (II -6)  and  using  Equations  (II-7)  through 

(II-9) 


'o  (‘M”') 

■ (7OR+l7OI) 


(JOR+1JOl)  + 


(11-10) 
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where 


and 


(I) 


2n 


(2n)  ! 


2n 

E 

j=i 


(ii-ii) 


Yoi=io  (“1)n 


(!) 


(2rtf-l ) 


(2nfl) ! 


(2n+l)  . ( 

A K 


(H-12) 


Introduction  of  the  second  boundary  condition  in  Equation  (11-10) 
requires  the  C2  = 0 because  YQ  (V^Vq  r)  approaches  infinity  at  r = 0. 

Making  use  of  the  first  boundary  condition  U(b,t)  = e1 Jt  or 
R(b,o>)  = 1 in  Equation  (II-4). 

1 = C1  J0  b)  = C1  (ber  *b  +ibeixb) 


(II-13) 


The  solution  of  Equation  (II-4)  that  is  the  frequency  response  function 
of  the  temperature  is  therefore 


R(r,u))  = 


ber  x + i bei 


ber 


+ i bei 


(11-14) 


Kelvin  functions  of  order  zero,  ber  x and  bei  x,  as  well  as  functions 
of  higher  order  are  tabulated  in  Reference  [5]  and  are  shown  in  Figure  1, 
The  tables  also  contain  Kelvin  functions  of  the  second  kind,  ker  x and 


kei  x which  are  related  to  YQ(x). 


The  higher  order  Kelvin  functions  ber  x and  bei  x needed  here 

m m 

are  also  expressed  in  series  form  [5] 


ber 


* -(p)  ilo 


and 


bei 


’■’fr) 


r , 1 

cos 

+ 

a 

k!  r (mfkfl) 


(H-15) 


sin 


k=0  k!  T (mfkfl) 


[(fa+  2O  Vi 


(H 


(11-16) 
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The  first  few  terms  of  ber^  x and  bei^  x are 

beri  x .JL,  (i  ) 

+ u2Y  i 

3!4!  \4  / * ' * 


(11-17) 


■‘,••4  [■  ■ i (s')  - rar  (s’)’ 


(11-18) 


The  magnitude  of  the  frequency  response  function  and  its  phase  angle 
may  be  obtained  from  Equation  (11-14)  as 


( ber  x ber  x,  + bei 
. v b 


x bei  x.  ) + (bei  x ber  x,  - ber  x bei  x, 

of ' b b > 


2 2 
ber  x.  + bei  x. 
b b 


bei  x ber  x,  - ber  x bei  x, 

tan  ?(r,j)  = r r — r~ r- : 

ber  x ber  x^  + bei  x bei  x^ 

Consequently  the  frequency  response  function  may  be  written  as 


(II-19) 


(11-20) 


R(r,j)  = j R(r , j) [ e 


(11-21) 


and  the  temperature  at  a point  as  a function  of  time  and  radial  position 
is  from  Equation  (II-2)  given  by 


U(r ,t)  = I R ( r , a) ) | e 


i(9+jt) 


The  real  part 


U(r,t)  = ! R ( r , o ) | cos  (9+  ot) 


(11-22) 


(11-23) 


indicates  the  time  delay  as  well  as  the  attenuation  of  the  temperature 
as  shown  in  Figure  2 for  a 22  in.  (55.9  cm)  propellant  cylinder  with 
steel  casing.  Figures  3 and  4 show  the  amplitude  and  the  phase  angle 
(time  delay)  in  the  interior  of  such  cylinders;  it  should  be  remembered 
that  both  |R(r,j)|  and  3(r,u)  are  functions  of  the  frequency,  u). 
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The  variance  of  the  temperature  at  a point  is  consequently  obtained 
from  Equations  (III-l)  and  (III-2)  by  substitution  of  Equation  (11-19) 
and  numerical  integration.  For  a narrow  band  random  excitation,  such 
as  that  produced  by  the  diurnal  cycle  ('j_  = 2itf_,  f_  = 1/24  hrs), 


Expanding  Equation  (11-19)  the  square  of  the  absolute  value  of  the 
frequency  response  function  simplifies  to 


If  the  input  temperature  is  considered  to  be  a band  limited  white 
noise  whose  one  sided  power  spectrum  is  a constant,  W , over  the  narrow 


frequency  range  from  f to  f„  and  is  zero  elsewhere  (Figure  5) 


where  frequency,  f,  in  cycles  per  hour  rather  than  radians  per  hour  is 
used. 


At  the  surface  the  variance  of  the  temperature  may  be  expressed  in 
terms  of  the  power  spectral  density.  Based  on  Equation  (III-2), 


I 

I 


1 


Hence  the  input  power  spectral  density  may  be  calculated  from  the 
surface  temperature  variance  as 


W 


\ 


E [Ui2(0)]  deg  / dee  C2  ) 
f 2 - f^  cycles/hr  \cycles/sec/ 


(HI-7) 


IV.  FREQUENCY  RESPONSE  FUNCTION  FOR  STRESS 

Thermal  stresses  in  long  cylinders  (plane  strain)  have  been 
developed  by  several  authors  such  as  Boley  and  Wiener  [7]  and  have  been 
calculated  for  a caseless  cylinder  subjected  to  sinusoidal  temperatures 
by  Dahl  [3],  A cylinder  with  a thin  case  has  been  treated  by  Williams, 
Blatz,  and  Schapery  [8]  for  arbitrary  temperatures.  Their  relations 
will  be  used  here  for  sinusoidal  and  random  temperature  distributions. 


' ‘ fe  f + (l-vKb2-a2,  (l*  i ) /"<r,rdr 

r a 

? / 


a 

(l-v 


U(r)rdr 


(IV-1) 


bf£ 


(i  + 4)  + — (i  + 4)  f 

-a  \ r / (l-v)  (b  -a  ) \ r 1 J 


aE  f 

+ y \ U(r)rdr 

(l-v)r  J 


a E U(r ) 
(l-v) 


= 2yb2p'  2 y aE 

az  ,.2  2,  ■r  w,2  2, 

(b  -a  ) (l-v)  (b  -a  ) 


U(r)rdr 


ayEU(r) 

(l-v) 


(IV-2) 


(IV-3) 


where  U(r)  is  the  difference  between  the  stress  free  (construction) 
temperature  and  its  momentary  value , a and  b are  the  inner  and  outer 

radii  of  the  cylinder  E,  v and  a are  the  modulus,  Poisson's  ratio  and 
a thermal  coefficient  of  expansion  for  the  material  and 
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E / U(r)rdr  - acU(b)  (l+vc>  E 

<b  ~a  > { 

(1+v)  b2  + + (l  - v 2) 


(IV-4) 


V,  -a 


E , vc>  and  a ! are  the  mechanical  parameters  of  the  case  materials  and 

h is  its  thickness, 
c 

For  a solid  cylinder  a = 0 and  Equations  (IV-1)  through  (IV-4)  are 
simplified  to 

S = - p1  + 7T7T  4 / rU(r)dr  - ~~  / rU(r)dr  (IV-5) 

( v)  [b  i r i 

~ 

S0  = - P'  + f rU(r)dr  + ~~  J rU(r)dr  - U(r) 

Lb  0 " 0 (IV-J6) 


S = 2V  pl  + I - 

z (1“v)  [b 


T-i  j rl)(r)rdr  - 

D < J 


( IV- 7 ) 


2q(1+ y) 


L b 


A* 


)rdr  - a£  (l+vc>  U(b) 


2 hF. 

(1+v)  (l-2v)  + (l-vc  ) 

c c 


. (IV-8) 


Gl.en  the  frequency  response  function  for  temperature,  R(r,oj;  of 
Equation  (11-14),  the  stress  response  functions  may  be  calculated  after 
the  integrals  in  Equations  (IV-5)  through  (IV-8)  are  evaluated 


J rU(r,oj)dr  = J r 


R(r,u)  dr 


(IV-9) 


Substituting  Equation  (11-14)  into  Equation  (IV-9) 


rU(r,<j)dr  = 


^ r [ber/|  r + i bei/|  r] 


her 

v rv 


dr  . (IV-10) 


- b + i bei/  - b 


The  integrals  of  Equation  (IV-10)  are  evaluated  individually. 
Changing  variables  to  x = r 


^ r ber  rdr  = ^ J °x  ber  x 
0 0 

j r bei  rdr  = ^ J ^x  bei  x 


(IV-11) 


(IV -12) 


From  Reference  [5] 


(1+m) 


1+  i / \ 

f dx  = - (f  - g . J 

m pr  \ nrt-1  nri-1  J 


(IV-13) 


where  m is  the  order  of  the  Kelvin  function  and  f and  g are  any  one  of 

m °m  J 

the  pairs  of  functions  given  below 


f = ber  x 
ro  m 


g = bei  x 
6m  m 


f = bei  x 
m m 


g = - ber  x 
m m 


(IV- 14) 


Hence  Equations  (IV-11)  and  (IV-12)  become  (m  = 0 is  not  usually  written: 
ber  qX  = ber  x) 

JT  berTf  rdr  = /|r  (be^  7|r  + be^  /§  r)  (IV-15) 


Jr  bei/|rdr=  /§r  (beij  /g  r + be^  /£  r ) . 


(IV-16 ) 
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At  r = 0 ber^  r = 0 and  bei^Y^  r = 0 
at  the  lower  limit  as  can  be  seen  from  the 
and  bei^x  [Equations  (11-15)  and  (11-16)]. 


at  r = 0,  the  express 


f 

0 


r ber/g 


(IV-5)  through  (IV-8)  approaches  the  value 
term  for  the  bei  function  approaches  zero. 
(IV-15)  and  (IV-16)  into  Equations  (IV-5), 
(IV-10)  and  using  the  abbreviations 


hence  both  integrals  vanish 
series  definitions  of  ber^x 
The  same  series  show  that 


r dr,  appearing  in  Equations 

of  +1/2  while  the  analogus 
Substituting  Equations 
(IV-6 ) , ( IV - 7 ) , (IV-8)  and 


(4^)  - 
(4')  * 


ber  J — r + i bei 


ys 

v a 


beri  r ■ beii^§ r 


(IV-17 ) 
(IV-18) 


(SI  r)  = berl  SI  r + beij  7f  r 


(IV- 19) 


2a  ( 1+ v ) / a 

lE  b 


[*(4>)  - AM . - 


P'(U))  = -• 


and 


x =Sz 


1ZM. 


+ ac  (l+vc), 


(1  + v)  (1  - 2v)  + (1  - v/)  ~~ 

c c 


(IV-20) 


Sr(r,j)  = 


- p'  + 


OE 

(1-v) 


< a 

2w 


b [k(xb^  4 i^^xb^  + r + **<«>]) 


h(xb) 


(IV-21) 


S (r,  .■)  - - p’  - 


(1-v) 


_£e_  fa_  )b  [k(V  + iix(b)J  + 2^¥h(x) 

(l-v)  * 2u>  | hvx^) 

(IV-23) 

where  Sr(r,u))  etc.  are  the  responses  to  a sinusoidal  surface  temperature 

of  unit  amplitude.  .It  has  been  assumed  in  the  previous  equations  that 
the  temperature  of  the  thin  casing  is  uniform  and  is  equal  to  the  ambient 
temperature. 


s2(r,u))  = 2v  ( p'  - 


At  the  center  of  the  cylinder  where  r = 0 the  stress  components  may 
be  obtained  from  Equations  (IV-5)  through  (IV-7)  by  utilizing  the 
limiting  values  of  the  integrals  discussed  previously.  Hence, 


where  p'  is  given  by  Equation  (IV-20). 


(IV-24) 


(IV-25) 


(TV-26) 


— 
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Equations  (IV-20)  through  (IV-23)  are  the  frequency  response 
functions  for  the  individual  stress  components.  They  are  presented  in 
Figures  6 and  7 for  the  material  properties  listed  in  Table  1.  The  time 
dependent  stresses  are  easily  obtained  by  multiplying  the  real  parts  of 

Equations  (IV-21)  through  (IV-23)  with  eiwt.  Figure  8 shows  the  time 
histories  of  and  at  various  positions  in  the  cylinder  as 


Si(r,u))  *»  |s^(r,co)|  cos  (ut  + <t>)i 


(IV-27) 


while  Figure  9 displays  the  variations  of  the  stress  components  through 
the  cylinder  at  fixed  times. 

V.  STRESSES  iN  THE  CASE 

It  has  been  assumed  in  the  foregoing  that  the  temperature  of 
the  case  is  uniform  through  its  thickness  and  is  equal  to  the  ambient 
temperature. 

The  case  stresses  are  given  in  Reference  [8]  for  the  general  plane 
strain  conditions  used  here  as 


i 2 / 2 \ aE  / u2\  r 

3r  = J¥LT'  f1  + T)+  C C 2 ' 2 f1  ' / rU<r>dr 

c (c  -b  ) ' r ' (1-v  ) (c  -b  ) \ r ' , 


a-vc)r  ' 


j rU(r)dr 


(V-l) 


,2  ' / 2\  aE  / ,2\  ,c 

, - -W  1 + -)  + s-f~r  i1*-)  l 

c (c  -b  ) ' r / (1-v  )(c-b2)  V r ’ C 

c 


rU(r)dr 


a E , a E u(r) 

c c I C C 7 

" ~2  ) rU<r>dr  - — ■ 

U-vc)r  b 


(V-2) 


2v  a E 


a v E U(r) 


2vb  D c c c f f V , c c c 

V - 7X727  + 7 XXTX  1 ■ XiXT 


z .2  2.  / N , 2 ,2.  J 

c (c  -b  ) (l-vc)(c  -b  ) 


a-vc) 


(V-3) 


where  c is  the  outer  radius  of  the  case.  As  a consequence  the  stress 
components  at  the  interface  between  case  and  cylinder  (r=b)  become: 


P' 


(&c 

\c2-b 


and 


(V-4) 


(V-5) 


(V-6 ) 


with  p'  given  by  Equation  (IV-4). 

Utilizing  the  material  properties  presented  in  Table  1 for  a 22  in. 
(55.9  cm)  propellant  cylinder  in  a 0.1  in.  (0.25  cm)  thick  steel  case 
the  interlaminar  stresses  in  the  case  are  evaluated  as: 


s 

r 

= 0.40  cos  ( 

'2fft 
k T 

0.92) 

1 (V-7 ) 

c 

se 

= 43.71  cos 

(¥ 

-0.92) 

/ 

1 (V-8) 

C 

S 

zc 

= 21.22  cos 

/ 2rtt 
V T 

-0.92) 

1 (V-9 ) 

where  T = 24  hrs  is  the  period  of  the  diurnal  cycle. 

VI.  MEAN  AND  VARIANCE  OF  THERMAL  STRESSES 

The  annual  thermal  cycle  consists  of  a yearly  mean  temperature 
on  which  a sinusoidal  variation  with  a frequency  of  one  per  year  is 
superimposed.  Thermal  stress  is  produced  by  the  difference  between  the 
stress  free  temperature  (the  temperature  at  the  time  of  construction) 
and  its  momentary  value.  If  it  is  assumed  that  the  cylinder  is  con- 
structed at  the  yearly  mean  temperature  then  only  the  slowly  varying 
annual  cycle  will  be  considered  in  the  determination  of  the  similarly 
varying  mean  values  of  the  stress  components.  While  these  could  be 
calculated  from  the  relations  developed  for  alternating  surface  tempera- 
tures, Equations  (IV-20)  through  (IV-23),  little  error  is  incurred  if 
it  is  assumed  that  the  temperature  throughout  the  cylinder  is  constant 
and  varies  sinusoidally  in  time  with  a frequency  of  1/8760  hrs.  It  can 
be  seen  from  Equations  (11-26)  and  (11-27)  and  the  series,  Equations 
(11-17)  and  (11-18),  that  the  temperature  at  the  center  of  a 22  in. 
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2 -6 

(55.9  cm)  diameter  propellant  cylinder  [a  = 1.27  in.  /hr  (0.228  x 10 
2 

m /s ) ] is  for  all  practical  purposes  unity  while  the  time  delay  is  of 
the  order  of  30  hrs , a small  amount  compared  to  the  yearly  period  of 
8760  hrs. 


Hence  the  annual  temperature  variation  throughout  the  cylinder  will 
be  characterized  as 


yr.t)  = uy  sin  (^-) 


(VI-1) 


Substituting  into  Equations  (IV-5)  through  (IV-8),  performing  the 
integrations  and  simplifying,  the  mean  values  of  the  stress  components 
are  obtained  as  follows: 


S = - p ' , S = - p ' and  S = 2vp  1 
rm  m 9m  m zm  m 


(VI-2) 


where 


E U a (1+v) 


a+v] 


m (l+v)(l-2v)  + (l-vc2)  ~§- 

c c 


f— ) 

\8760  / 


(Vl-3) 


Equations  (V-2)  and  (V-3)  may  also  be  derived  as  a limiting  case  of  the 
sinusoidal  thermal  input  by  permitting  the  frequency,  to,  to  assume  very 

small  values  ^co  = f°r  the  annual  cycle^  . Utilizing  the  first 

two  terms  of  the  series  expansions  of  ber^,  bei^,  ber  and  bei  [Equations 

( 1 1 - 7 ) , (II-8),  (."1-17),  and  (11-18)]  the  integrals  of  Equations  (IV-15) 
and  (IV- 16)  reduce  to 


r ber 


rz  . r“ 
N/-  r dr  = t 


(VI-4) 


I r bei'/f  r dr  = a 16 


(VI-5  ) 


and  the  function  h of  Equation  (IV-17)  when  higher  order  terms  containing 
to  are  dropped  becomes 
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(VI-6) 


= 1 + 


2 

a)  r_ 

a 4 


i 


Substituting  these  expressions  into  the  equations  for  the  stresses  and 
neglecting  the  imaginary  terms  (they  alone  contain  w)  Equations  (VI-2) 
and  (VI-3)  are  again  obtained. 

The  value  of  p^'  has  been  calculated  for  the  material  properties 
shown  in  Table  1 as 

'.'"•‘’‘Vlw)  • (VI-7) 


If  there  is  a difference,  Um>  between  the  stress  free  temperature  and 

the  yearly  mean  an  additional  term  analogous  to  Equation  (VI-3),  but 

with  U replacing  U sir.  ( — \ , ) would  be  added  to  Equation  (VI-3), 
m y \o/ou/ 


The  narrow  band  random  stresses  produced  by  the  diurnal  cycle  are 
superimposed  on  these  slowly  varying  mean  values.  The  variances  of  the 
stress  components  may  be  calculated  in  a manner  analogously  to  the 
temperature  variance  as  shown  in  Section  III. 


The  power  spectrum  of  any  one  of  the  stress  components  is  given 
directly  from  the  power  spectrum  of  surface  temperatures.  From 
Equation  (III-l) 


Wsj(f)  = | Sj (r, f ) | 2 W (f) 


j = r,0,z 


(VI-8) 


While  the  variance  of  a stress  component  may  be  obtained  by  numerical 
integration  as 

00 

E [sj2(r)]  = / I Sj  (r » f ) 1 2 Wu  < f > df  • (VI-9) 

0 1 


For  the  band  limited  white  noise  surface  temperature  input  with  centnl 
frequency  expressed  in  cycles  per  hour  and  with  (f)  = constant,  the 

2 2 

one  sided  power  spectral  density  in  (°F)  /cy/hr  or  (°C)  /cy/hr  is 


(VI-10) 
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The  standard  deviations  which  are  the  square  roots  of  the  variances 
have  been  calculated  with  a rectangular  band  limited  white  noise  centered 
over  the  24-hour  period  as  shown  in  Figure  5. 


In  this  narrow  range  the  frequency  response  functions  are  relatively 
constant.  As  a consequence  the  standard  deviations  of  the  temperature 
and  the  stress  components  for  a 1°F  (0.555°C)  surface  temperature  stand- 
ard deviation  follow  the  curves  of  Figure  3 within  plotting  accuracy. 

In  the  particular  geometry  and  material  properties  chosen,  the  curves 
for  the  amplitudes  of  the  frequency  response  functions  also  indicate 
the  standard  deviations  of  the  various  functions  for  a one  degree 
Fahrenheit  surface  temperature  standard  deviation. 


The  average  frequency  of  the  response  to  a narrow  band  excitation 
may  be  obtained  from  random  process  theory  [6], 


According  to  Equation  (VI-8)  the  average  frequency  of  the  stress  compo- 
nents at  any  point  within  the  cylinder  may  then  be  expressed  in  terms 
of  the  frequency  response  function  and  the  input  power  spectral  density 
as 


When  the  input  power  spectral  density  is  a constant 


df 


J2  f2  |s.(r,f)|2 


f 2 = -i-r . (VI-13) 

f 

/ | Sj (r ,f ) | 2 df 

fl 

Integration  of  Equation  (VI-12)  for  the  narrow  frequency  range  of 
f^  = 1/23  cy/hr  and  f^  = 1/25  cy/hr  shows  that  the  average  frequency 


of  the  output  stress  components  is  very  close  to  f = 1/24  because  the 
frequency  response  functions  are  nearly  constant  in  this  range. 


i VII.  DISTRIBUTION  OF  TEMPERATURE  PEAKS 


If  the  ambient  surface  temperature  in  addition  to  being  a 
narrow  band  random  process  has  Gaussian  probability  of  exceedance  the 
distribution  of  amplitude  peaks  may  be  calculated.  The  density  function 
of  amplitudes  follows  a Rayleigh  distribution  [6]  of  the  form 


(VII-1) 


where  a is  the  random  amplitude  of  the  diurnal  cycle  and  is  the 
standard  deviation  of  hourly  temperature  measurements. 

The  probability  of  encountering  an  amplitude  greater  than  a given 
value  a is  the  integral  of  Equation  (VII-1) 


P[A  > a] 


(VII-2  ) 


The  mean  amplitude  of  the  process  [9]  is 

I -M  a - 1.25  a (VII-3) 

2 u u 

. 
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while  the  standard  deviation  of  amplitudes,  a , may  be  given  in  terms  of 

a 

the  standard  deviation  of  hourly  temperature  readings  cr^  as 


0.65  a 


(VII-4) 


Expressions  analogous  to  Equations  (VII-1),  (VII-2),  and  (VII-3)  may  be 
written  for  the  amplitude  distributions  of  stresses  and  temperatures  at 
various  positions  in  the  cylinder  by  substituting  the  appropriate 
deviations,  a , a , a , and  o in  these  equations. 

Sr  S9  S z U 


The  values  of  the  required  standard  deviations  for  a one  degree 
surface  temperature  standard  deviation  may  be  obtained  from  Figure  3. 

VIII.  SAMPLE  CALCULATIONS 

In  order  to  calculate  the  stress  response  of  a structure  to 
thermal  variations,  temperature  records  are  necessary.  The  information 
required  consists  of  the  average  temperature  at  the  location  of  the 
structure,  the  amplitude  of  the  annual  cycle,  and  the  standard  deviation 
of  hourly  temperatures. 

Assuming  that  the  solid  propellant  was  cast  at  160°F  (71.1°C)  and 
is  deployed  at  a location  where  the  annual  mean  temperature  equal 
60°F  (15.6°C),  the  amplitude  of  the  annual  cycle  (U^)  = 25°F  (13.9°C) 

while  the  standard  deviation  of  hourly  temperature  measurements 
(a  ) = 8°F  (4.4°C) . 

Using  Equations  (VI-7),  (VII-3),  and  Figure  6,  the  tangential 
stress  at  the  interface  between  propellant  and  steel  case  will  consist 
of  the  following  components;  the  slowly  varying  mean  stress 

S,  = -1.171  [25  sin  - (160  - 60)] 

and  the  average  amplitude  of  the  diurnal  cycle 

S (b)  = 1.25  a x a = 1.25  X 8 x 0.41  = 4.1  psi  (28.3  KPa) 


whose  standard  deviation  is 


a,  (b)  = 0.65  o X a = 2.13  psi  (14.7  KPa) 
9 us,  r 

a 9 


Consequently  at  summer  temperatures  the  interfacial  tangential  stress 
will,  on  the  average,  cycle  between  83.7  (577)  and  121.2  psi  (836  KPa) 
while  during  the  winter,  stresses  between  113  (779)  and  150.5  psi 
(1038  KPa)  will  be  observed. 


With  ±3  standard  deviations  these  stresses  will  have  maximum  and 
minimum  values  of  77.3  (533)  and  156.9  psi  (1082  KPa). 

At  an  interior  point  stresses  may  be  calculated  in  a similar  manner 
if  the  appropriate  standard  deviations  are  used  from  Figure  6. 

IX.  CONCLUSIONS 

Thermal  stresses  produced  by  a narrow  band  random  surface 
temperature  variation  in  a solid  propellant  cylinder  encased  in  a steel 
shell  have  been  calculated. 

A computer  program  has  been  developed  for  the  calculation  of  stress 
response  to  sinusoidal  and  random  thermal  inputs. 

It  has  been  found  that  the  diurnal  temperature  cycle  produces  cyclic 
stresses  whose  amplitude  remains  essentially  constant  throughout  the 
cylinder.  This  is  in  contrast  to  the  variations  of  the  temperature 
which  are  attenuated  at  the  interior. 

It  is  expected  that  this  behavior  is  geometry  dependent  and  that 
motors  with  smaller  diameters  would  exhibit  different  characteristics. 

For  the  large  diameter  motor  investigated  here  the  effects  of  the 
diurnal  cycle  are  insignificant  compared  to  those  produced  by  the  annual 
cycle. 

In  a subsequent  study  stress-strength  interference  and  probability 
of  failure  of  such  motors  may  be  examined. 
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Figure  5.  Band  limited  white 
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Figure  6.  Amplitudes  of  frequency  response  functions,  standard 
deviations  of  response  functions  for  unit  surface  temperature 
[b  = 11  in.  (27.9  cm)]. 
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